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摘 要 三 次 多 项 式 温度 分 布 模型 反映 了 液 滴 内 部 基于 希 尔 球 涡 的 实际 流动 ， 可 以 较为 准确 地 预测 单个 液 滴 蒸 发 . 鉴于 
三 次 多 项 式 温度 分 布 模型 在 预测 单个 液 滴 蒸发 的 准确 高 效 性 ,本文 以 n-Decane 油 滴 蒸发 为 例 , 采用 三 次 多 项 式 模型 研究 
了 通过 改变 初 温和 初始 气流 速度 来 改变 初始 雷诺 数 对 液 滴 蒸 发 的 影响 , 结果 表明 通过 改变 气流 速度 来 改变 雷诺 数 的 情况 ， 
不 同 雷 诺 数 下 液 滴 燕 发 率 与 内 部 温度 分 布 存 在 较 大 差异 ， 以 n-Decane, n-Heptane 和 n-Dodecane 三 种 不 同 黏度 的 油 滴 
燕 发 为 例 , 对 比 了 液 滴 条 度 对 蒸发 的 影响 , 发 现 不 同 黏度 液 滴 的 蒸发 率 与 内 部 温度 分 布 差异 较 大 , 条 度 越 小 , 燕 发 率 越 快 ; 
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Abstract A 3rd_order Wo 
proposed recently, in Wi the enhancement of heat transport by A circulation was correctly 


omial temperature profile liquid model for droplet evaporation was 


reflected due to the introduction of Hill’s vortex. Thus it could pfeghict droplet temperature and 


evaporation rate well. Based on this well-built liquid evaporati bdel, this paper investigated the 


effect of initial parameters including initial droplet temperatuie and ambient gas velocity on droplet 


evaporation rate. It was found that the initial gas vio a larger effect on droplet evaporation 


rate and temperature distribution than the initial NA emperature. Further, several kinds of liquid 
Ne 


droplets (n-Decane, n-Heptane and n-Dodecane employed to study liquid viscosity on droplet 


evaporation. The result indicated that a smaller. scosity of liquid droplet could lead to a larger droplet 
evaporation rate and a greater temperatur& tadient inside droplet in the transient evaporation period. 
Key words droplet evaporation; hill’s vortex; 3rd-order polynomial temperature profile model; liq- 


uid viscosity 
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0 引言 


液 滴 蒸发 涉及 液 滴 与 周围 环境 的 传 热 、 传 质 以 
及 动量 交换 ， 是 影响 整个 喷雾 液 滴 群 喷雾 动力 学 与 
热 特性 行为 的 重要 因素 ,同时 在 工业 应 用 中 也 非常 
重要 . 例如 , 液 滴 蒸发 对 内 燃 机 燃油 的 燃烧 效率 , 激 
光 皮 肤 手 术 中 制冷 剂 喷雾 冷却 在 皮肤 表面 的 冷却 效 


有 重要 理论 及 实际 意义 , 但 是 既 要 保证 模型 的 准确 
性 又 要 保证 计算 的 高 效 性 却 往 往 是 一 个 相互 矛盾 的 
问题 ,最 近 本 课题 组 提出 一 种 基于 液 滴 内 部 温度 服 
从 三 次 多 项 式 分 布 的 液 相模 型 ， 考虑 了 基于 希 尔 球 
涡 的 实际 流动 对 液 滴 内 部 传 热 与 蒸发 的 影响 ， 即 当 


率 以 及 冷却 塔 的 换 热 效果 等 均 具 有 重要 影响 Pl. 
因此 ， 建 立 准确 并 且 计 算 高 效 的 液 滴 蒸发 模型 具 
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气流 流 过 液 滴 表面 ， 作 用 于 液 滴 表 面 的 前 切 力 会 使 
液 滴 产生 一 个 内 部 流动 ， 从 而 加 强 了 液 滴 内 部 的 传 
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热 回 . 与 集 总 参数 模型 和 二 次 多 项 式 温度 分 布 模型 
相 比 ， 新 的 三 次 多 项 式 模型 可 以 更 为 准确 地 反映 液 
滴 内 部 的 实际 流动 情况 与 温度 梯度 "站 ， 与 导热 模 
型 和 有 效 导 热 模 型 相 比 ， 三 次 模型 预测 结果 也 更 为 
准确 ， 且 计算 简单 ， 效 率 更 高 由。 三 次 多 项 式 温度 
分 布 模型 平衡 了 针对 单个 液 滴 蒸 发 的 液 相 模型 的 准 
确 性 与 高 效 性 ， 适 合 被 广泛 地 应 用 于 喷雾 CFD 计 
算 中 . 

鉴于 三 次 多 项 式 温度 分 布 模型 在 预测 液 滴 蒸发 
过 程 中 的 优越 性 ， 本 文 将 采用 三 次 多 项 式 温度 分 布 
模型 对 单个 液 滴 莹 发 过 程 进行 研究 . 以 n-Decane 油 
滴 蒸 发 为 例 ， 采 用 三 次 多 项 式 模型 研究 初始 条 件 对 
液 滴 蒸发 的 影响 ， 分 析 通 过 改变 初 温和 初始 气流 速 
度 来 改变 初始 雷诺 数 对 液 滴 蒸 发 的 影响 。 由 于 液 滴 
条 度 对 其 流动 具有 重要 影响 ， 本 文 以 n-Decane, n- 
Heptane 和 "-Dodecane 三 种 不 同 务 度 的 油 滴 营 发 为 
例 ， 应 用 新 模型 对 比 研 究 了 液 滴 黏 度 对 蒸发 及 内 部 
传 热 的 影响 ， 


1 省 滴 莹 发 液 相模 型 

如 果 不 考虑 液 消 内 部 温度 梯度 ， 把 深 注 滤 休 当 
作为 一 个 温度 处 理 ， 即 集 总 参数 模 现 消 能 时 
衡 方程 为 : 0 ~ 


1 
h(Ts 一 Ts) = 了 rsaiee Er piLrs (1) 
式 中 , h, Ts 与 工分 别 表示 对 流 换 热 系数 ,环境 气 
体温 度 和 液 滴 表面 温度 ， Ts 为 液 滴 半径 ， DO，Cp,1， L 


和 六 分 别 代表 液体 密度 ,液体 比热容 ， 节 发 潜 热 和 
液 滴 半径 变化 率 . 


虽然 集 总 参数 模型 计算 效率 高 ， 但 是 不 能 反映 、\ 
” 应 的 环境 气体 速度 分 别 为 0 m/s, 1 m/s, 3.8 m/s 和 


液 滴 内 部 的 温度 分 布 .而 近来 由 本 课题 组 提出 此 
于 三 次 多 项 式 温 度 分 布 和 液 滴 内 部 基于 希 尔 球 滴 的 
实际 流动 的 液 相 模型 不 仅 可 以 考虑 液 滴 内 部 温度 分 
布 69- 韦 ， 并 且 计 算 简单 高 效 . 三 次 模型 的 基本 思 
路 是 用 一 个 三 次 多 项 式 函 数 来 近似 表示 液 滴 内 部 温 
度 分 布 ， 用 平均 温度 来 代替 集 总 参数 模型 当中 的 液 
滴 整 体温 度 , 由 平均 温度 了 表示 能 量 方程 为 : 


dT 二 二 L h 
tb 
S TsplCp,l 


— 


根据 液 滴 表面 能 量 平 衡 、 液 滴 中 心 温度 对 称 、 希 
尔 涡 沾 止 条 件 以 及 平均 温度 定义 式 四 项 条 件 ， 可 
以 求 得 三 次 多 项 式 温度 分 布 模型 的 表面 温度 到 与 
平均 温度 的 关系 式 可 以 表示 为 (详细 推导 可 参考 
文献 [6]): 


E> 
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_ T+0.043XMgNuT, 
?1+0.043MNu 


piLrsTs 
和 (1+0.043Mg1Nu) 


式 中 , 和 gj 为 液 滴 表面 混合 气体 导热 系数 与 液体 导热 
系数 之 比 . 式 (2)、(3) 给 出 了 液 滴 蒸发 过 程 中 基于 希 
尔 球 涡 理 论 的 三 次 多 项 式 温度 分 布 液 相模 型 的 完整 
表达 ， 

此 外 ,关于 决定 液 滴 表 面 与 环境 气体 间 传 热 传 
质 的 Nu 和 Sh 气相 模型 , 本 文中 选用 文献 [12] 提出 
的 适合 于 不 同 蒸发 率 下 N-G-R-M 气相 模型 . 


2 三 次 模型 计算 液 滴 节 发 结 来 与 讨论 
2 ale 
次 多 项 式 温度 分 布 模型 在 预测 液 滴 蒸 发 


0 
滴 


0.086 (3) 


加 热 蒸发 的 影响 , Wong 与 Lin 的 实验 结果 表 
月 ， 在 雷诺 数 0~100 的 范围 内 , 雷诺 数 的 改变 对 液 
滴 燕 发 的 影响 很 小 ,但 是 只 考虑 了 改变 初始 温度 来 
改变 雷诺 数 的 情况 "3 引 . 而 实际 上 雷诺 数 随 液 滴 初 始 
直径 、 周 围 环境 气体 流速 和 液 滴 表 面 混 合 气 体 的 热 
物性 的 变化 而 改变 ,并且 液 滴 内 部 流动 与 传 热 受 周 
围 气 流速 度 的 影响 更 大 。 因 此 为 了 进一步 研究 初始 
雷诺 数 对 液 滴 蒸发 的 影响 ， 本 文 分 别 考虑 了 通过 改 


变 初 始 温 度 和 周转 环 境 气 体 流 速 来 改变 初始 雷诺 数 
的 情况 ， 并 允 湾 滴 蒸 发 率 和 内 部 温度 分 布 随时 间 变 


化 进行 陛 较 . 在 2.1.1 中 , 固定 环境 气体 流速 为 1 
m/sx 狂 诺 数 为 17，60 和 100 时 所 对 应 的 初 温 分 别 
类 入 8K, 307K 和 301K 在 2.1.2 中 ,固定 液 滴 初 
温度 为 315K, 雷诺 数 为 0，17，60 和 100 时 所 对 


6.3 m/s。 在 两 种 情况 下 ， 液 滴 初始 直径 都 为 2 mm， 
周围 气体 温度 为 1000 K. 
2.1.1 通过 改变 初始 温度 改变 初始 雷诺 数 

图 1 为 三 次 模型 计算 w-Decane 油 滴 在 不 同 初始 
温度 下 液 滴 直 径 随时 间 的 变化 . 当 液 滴 初始 温度 较 
高 时 ，"-Decane 液 滴 的 蒸发 率 也 较 高 ， 但 在 不 同 初 
始 雷 诺 数 下 计算 的 液 滴 蒸 发 率 没有 较 大 的 差异 . 初 
始 温度 的 增加 导致 雷诺 数 Re 与 施 密 特 数 Sc 的 减 
小 ， 从 而 导致 舍 伍 德 数 5h 的 减 小 . 同时 ,扩散 系数 
减 小 . 因此 , 在 雷诺 数 Res<O(10)~O(100) 的 范围 内 ， 
液 滴 燕 发 率 没有 明显 的 差异 ,与 Wong 和 Lin 之 前 
的 结论 一 致 . 图 2 为 三 次 模型 计算 n-Decane 油 滴 在 
不 同时 刻 不 同 初始 温度 下 内 部 温度 分 布 情况 . 液 滴 
蒸发 的 初始 阶段 , 在 雷诺 数 Re=<O(10)~O(100) 的 范 
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围 内 ,不 同 雷诺 数 下 计算 得 到 的 液 滴 内 部 温度 分 布 


有 一 定 的 差异 , 但 是 随 着 时 间 增 加 , 差异 逐渐 减 小 . 


v=] m/s 

Re=17, 70=313 K 
一 -Re=60, 70=307 K 
一 - Re=100,70=301.6 K 


0 2 4 6 
t/s 
图 1 三 次 模型 计算 n-Decane 油 滴 在 不 同 初始 温度 下 液 滴 
直径 随时 间 的 变化 
Fig. 1 Variation of n-Decane droplet diameter versus time at 


different initial temperature 


图 2 
时 刻 (t=0.4 s 和 t=2.4 s) 内 部 温度 分 布 情况 


Fig. 2 JInternal temperature distribution of n-Decane droplet 


with varied initial temperature at different time 
(t=0.4 s and t=2.4 s) 


2.1.2 通过 改变 气流 速度 改变 初始 雷诺 数 

图 3 为 不 同 气流 速度 时 三 次 模型 计算 n-Decane 
油 滴 的 直径 随时 间 变 化 . 可 以 看 到 , 通过 改变 气流 速 
度 来 改变 初始 雷诺 数 的 情况 与 通过 改变 初始 温度 来 
改变 初始 雷诺 数 的 情况 有 较 大 区 别 . 不 同 初始 雷诺 
数 下 液 滴 蒸发 率 存 在 较 大 差异 . 初始 温度 固定 , 气流 
速度 越 大 , 初始 雷诺 数 越 大 , 液 滴 蒸 发 越 快 . 根据 希 
尔 球 涡 理论 , 流 过 液 滴 表面 的 气体 流速 越 大 , 作用 于 
液 滴 表面 的 剪 切 力 所 产 生 的 内 部 流动 得 到 加 强 ， 因 
而 促进 了 液 滴 蒸发 . 当 气 流速 度 为 零 时 (v = 0 m/s)， 
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即 液 滴 与 环境 气体 相对 静止 时 ， 液 滴 蒸发 很 慢 。 可 
以 得 出 ， 液 滴 内 部 流动 对 液 滴 加 热 与 殖 发 有 很 大 的 
影响 , 而 液 滴 内 部 流动 受 环境 气流 速度 影响 较 大 . 


T=315 K 
Re=0,v=0 m/s 
一 -Re=17,v=1 m/s 


一 - Re=60,v=3.8 m/s 
--- Re=100,v=6.3 m/s 


ts 
AAA 
CS n-Decane 油 滴 的 直径 
随时 间 变 化 


Variation of n-Decane droplet diameter versus time at 


各 


图 4 为 三 次 模型 计算 n-Decane 油 滴 在 不 同 气 流 
速度 下 不 同时 刻 内 部 温度 分 布 情况 .在 液 滴 蒸发 的 
初始 阶段 (t = 0.4 s), 当初 始 雷诺 数 为 17, 60 和 100 
的 情况 下 (v = 1 m/s, 3.8 m/s, 6.3 m/s), n-Decane 
液 滴 内 部 具有 较 大 的 温度 梯度 .由 于 外 界 气流 作用 
于 液 滴 表 面 的 剪 胡 为 所 产生 的 内 部 流动 加 强 了 液 泣 
表面 向 中 心急 传 热 , 因此 气流 速度 越 大 , 液 滴 中 心 的 
加 热 越 明 最 表面 温度 与 中 心 温度 都 明显 高 于 他 们 

Ho 因此 液 滴 内 部 的 温度 梯度 越 大 . 随 
0 内 部 温度 梯度 减 小 , 气流 速度 越 大 ， 
温度 越 早 趋 于 均匀 分 布 . 可 以 看 出 , 在 t=2.4 s 


different ambient gas velocity 


三 次 模型 计算 n-Decane Cs 时 ,对 于 雷诺 数 为 60, 100(v= 3.8 m/s, 6.3 m/s) 的 


情况 , 液 滴 达 到 稳 态 蒸发 阶段 . 而 对 于 雷诺 数 为 0(v 
= 0 m/s) 的 情况 , 液 滴 整 体温 度 较 低 ， 并 且 在 液 滴 
蒸发 的 整个 过 程 中 , 内 部 温度 总 是 趋 于 均匀 分 布 , 这 
是 由 于 当 液 滴 与 环境 气体 相对 静止 时 , 液 滴 不 存在 
基于 希 尔 球 涡 的 内 部 流动 ， 
2.2 液 滴 黏 度 对 蒸发 的 影响 

基于 希 尔 球 涡 理 论 的 三 次 多 项 式 温度 分 布 模 
型 考虑 了 液 滴 内 部 温度 分 布 梯度 ， 反 映 了 液 滴 内 
部 的 实际 流动 情况 ， 当 气流 流 过 液 滴 表 面 ， 作 用 
于 液 滴 表面 的 前 切 力 会 使 液 滴 产生 一 个 内 部 流动 ， 
从 而 加 强 了 液 滴 内 部 的 传 热 过 程 。 对 于 不 同 种 类 
的 油 滴 ， 灰 度 的 不 同 会 对 液 滴 内 部 流动 情况 产生 
影响 ， 因 此 本 节 采 用 三 次 多 项 式 模 型 并 选取 了 黏 
度 不 同 的 三 种 油 滴 (nmn-Heptane，n-Decane 和 7- 
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Dodecane) 来 研究 液 滴 务 度 对 闵 发 的 影响 。 其 中 rm 
heptane 液 滴 的 黏度 为 7.56x10-4uPa.s, 沸点 温度 为 
98.4°C; n-Decane 液 滴 的 黏度 为 1.28x10-3 Pa.s， 
沸点 温度 为 174*C;，n-Dodecane 液 滴 的 条 度 为 
4.22x10-2 uPa.s， 沸 点 温度 为 215.9*C。 计算 中 ， 
周围 气体 温度 为 741 K， 压 强 为 0.1 MPa， 液 滴 初 
始 温 度 为 300 K, 液 滴 初 始 直 径 为 0.6 mm。 


Re=0 Re=17 Re=60 Re=100 


. 全 v 


T=315 K 


0.0 0.2 0.4 
Kr. 


图 4 三 次 模型 计算 n-Decane 油 滴 在 不 同 气 流速 度 及 2 
(t=0.4 s, t=2.4 s) a 


Fig. 4 Internal temperature distribution of NY VY 


exposed to varied gas velocity at times Nh “s and t=2.4 s) 


党 1 
图 5 为 雷诺 数 为 17 时 六 次 型 计算 n-Heptane,， 
n-Decane 和 Dodecane 油 滴 的 蒸发 率 , 内 部 温度 分 
布 和 表面 温度 随时 间 变 化 . 与 其 他 两 种 液 滴 相 比 , n- 
Dodecane 液 滴 的 蒸发 率 最 低 , 这 是 由 于 n-Dodecane 
较 高 的 液 滴 称 度 造 成 的 . 液 滴 务 度 越 大 , 对 内 部 流动 
的 阻碍 就 越 大 , 从 而 减弱 了 液 滴 热量 传递 , 导致 蒸发 
率 偏 小 . 当 雷 诺 数 为 17 时 , 在 1/d2 = 0.1 的 时 刻 ， 可 
以 明显 看 到 三 种 液 滴 内 部 存在 温度 梯度 ， 并且 
高 备 度 液 滴 阻碍 内 部 涡 旋 流动 ,所 以 当 液 滴 刍 度 越 
小 , 液 滴 内 部 流动 越 强 , 对 中 心 温度 的 加 热 促进 作用 
越 明 显 ， 因 此 内 部 温度 梯度 越 大 . 在 2/d = 1 的 时 
刻 , 由 于 n-Heptane 油 滴 已 达到 稳 态 蒸发 阶段 , 内 部 
温度 分 布 趋 于 均匀 ,. 而 n-Decane 油 滴 和 n-Dodecane 
油 滴 内 部 仍然 存在 一 定 的 温度 梯度 。 从 图 5(c) 中 可 
知 ， 在 液 滴 蒸 发 的 初始 阶段 ，n-Heptane 油 滴 由 于 
笑 度 最 小 ,表面 温度 上 升 最 快 ， n-Decane 油 滴 与 nm- 
Dodecane 油 滴 的 表面 温度 上 升 速率 相近 , n-Heptane 
油 滴 大 约 0.2 s 达到 稳 态 温度 ，n-Decane 油 滴 大 约 
1.5 s 达到 稳 态 温度 ， 而 mn-Dodecane 达到 稳 态 温度 
所 需 时 间 最 久 ， 大约 2.5 s。 三 种 油 滴 最 后 达到 的 稳 
态 温度 的 大 小 关系 与 它们 的 沸点 温度 的 大 小 关系 相 
一 致 ， 
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(b) 液 滴 内 部 温度 


(b) Internal temperature distribution 
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图 5 三 次 模型 计算 n-Heptane, n-Decane 和 n-Dodecane 油 滴 
的 燕 发 率 、 内 部 温度 分 布 和 表面 温度 随时 间 变 化 (Reo=17) 
Fig. 5 Evaporation characteristics of n-Heptane, n-Decane 

and n-Dodecane droplets (Reo=17) 


3 结 论 


本 文采 用 基于 希 尔 球 涡 理 论 的 三 次 多 项 式 温度 
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分 布 液 相 模型 针对 喷雾 过 程 中 的 单个 液 滴 蒸发 过 程 
开展 理论 研究 . 首先 以 mn-Decane 油 滴 蒸 发 为 例 ,分 
析 了 分 别 通过 改变 液 滴 初 始 温度 与 外 界 气流 速度 来 
改变 初始 雷 详 数 对 液 滴 加 热 与 燕 发 的 影响 ， 结 果 表 
明 ， 当 通过 改变 初始 温度 来 改变 初始 雷诺 数 时 ， 在 
不 同 雷诺 数 下 液 滴 蒸 发 率 情况 相近 ， 

而 通过 改变 气流 速度 来 改变 初始 雷诺 数 对 液 滴 蒸发 
率 与 温度 分 布 有 显著 的 影响 ， 气流 速度 的 增加 引起 
内 部 流动 的 加 强 , 导致 液 滴 蒸发 变 快 , 在 蒸发 的 初始 
阶段 内 部 温度 梯度 较为 明显 。 接 着 以 n-Decane,，n- 
Heptane 和 mn-Dodecane 三 种 不 同 黏度 的 油 滴 蒸 发 为 
例 ， 对 比 了 不 同 竺 度 液 滴 的 蒸发 率 与 温度 分 布 的 差 
异 , 液 滴 黏 度 的 增 大 导致 了 内 部 流动 和 传 热 的 减弱 ， 

因此 液 滴 黏 度 越 大 ， 蒸 发 越 慢 ， 在 液 滴 蒸 发 的 初始 
阶段 内 部 温度 分 布 梯度 越 不 明显 . 因此， 三 次 多 项 
式 温 度 分 布 模型 在 对 单个 液 滴 蒸 发 的 计算 中 体现 出 
了 实际 的 内 部 流动 对 液 滴 蔡 发 与 加 热 的 促进 作用 . 
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